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We discuss the localization of wavefunctions along planes 
containing the shortest periodic orbits in a three-dimensional 
billiard system with axial symmetry. This model mimicks 
the self-consistent mean field of a heavy nucleus at deforma- 
tions that occur characteristically during the fission process 
mM ■ Many actinide nuclei become unstable against left-right 
asymmetric deformations, which results in asymmetric frag- 
ment mass distributions. Recently we have shown |3|y that 
the onset of this asymmetry can be explained in the semi- 
classical periodic orbit theory by a few short periodic orbits 
lying in planes perpendicular to the symmetry axis. Presently 
we show that these orbits are surrounded by small islands of 
stability in an otherwise chaotic phase space, and that the 
wavefunctions of the diabatic quantum states that are most 
sensitive to the left-right asymmetry have their extrema in 
the same planes. An EBK quantization of the classical mo- 
tion near these planes reproduces the exact eigenenergies of 
the diabatic quantum states surprisingly well. 

PACS numbers: 24.75.+i, 03.65. Sq, 21.10.Dr, 47.20. Ky 



We have recently applied SM] the periodic orbit the- 
ory [pP] to a three-dimensional cavity model with axially 
symmetric deformations that typically occur near the iso- 
mer minimum and the second maximum of the character- 
istic double- humped fission barrier ||l|,y| of many actinide 
nuclei. The boundary of the cavity in cylindrical coor- 
dinates {z,p) is given by a shape function p{z); the de- 
formations are defined by an elongation parameter c and 
an octupole-type left-right asymmetry parameter a (see 
Ref. ^ for the detailed definitions; the neck parameter 
h is fixed to be zero in the present work). The shortest 
periodic orbits are found in planes perpendicular to the 
symmetry (z) axis at locations Zi given B by p'{zi) — 0; 
they are just the shortest polygons (diameter, triangle, 
square, etc.) inscribed into the circular cross sections 
of the billiard system with radii p{zi). The oscillating 
part SE of the total energy of the system containing N 
fermions (we do not distinguish neutrons from protons 
and neglect the Coulomb and spin-orbit interactions) was 
calculated by the appropriate semiclassical trace formula. 
A uniform approximation was introduced to describe the 
bifurcation of a single orbit plane (at zq = 0) into three 
orbit planes Zi {i — 0, 1,2) at the moment of the neck 
formation where the shape function p{z) becomes convex 
(see Ref. [pi for details). The gross-shell features in SE 
were emphasized by convolution of the trace formula over 



the wave number khy a, Gaussian of width 0.6/i? (i? is the 
radius of the spherical nucleus), producing a smearing of 
the shell structure similar to that caused by the residual 
pairing interaction in realistic nuclear models p| . 

In the quantum-mechanical calculations for the energy 
shell correction SE, derived from the quantum spectrum 
of realistic nuclear shell-model potentials using Strutin- 
sky's shell-correction method 0, it was found that the 
outer fission barrier of many actinide nuclei is unstable 
against octupole-type deformations ||8||. In our present 
parametrization, this barrier is located at c = 1.53; by 
varying a from zero to ^ 0.13, the energy SE is lowered 
by about ~ 2.5 — 3 MeV. (For the present qualitative dis- 
cussion, we neglect the smooth liquid-drop model part of 
the total energy which varies much less in this restricted 
region of deformations.) This result could be well repro- 
duced in our semiclassical calculations MM- Thereby, it 



/>(z) 



0.5 
0.0 
-0.5 





-1.5 



-1.0 -0.5 0.0 0.5 1.0 1.5 
z 



1.0 



0.5 



0.0 



-1.0 



■# fr 


' ■;•: A . w ^- ' ^ 






'■ - 1 


<■>.>. V 








r" ^ i 






> ^ 'f' 




y4 e,i- - 








n. Si 


fx---^r.i 




'i^j* 


.,-'.'.? 


;'^ i 








1 ^ 


•«' l-fV''.' 


^^'^ 


^^-^ 




' -^ 









w 



FIG. 1. Top: Cavity shape p{z) at elongation c — 1.53 and 
asymmetry a — 0.0. The vertical lines indicate the planes 
containing the shortest periodic orbits. Bottom: Poincare 
surface of section (p, (fi) for L^ = 0. See the text for details. 



turned out to be sufRcient to include the two shortest 
periodic orbits (diameters and triangles) in each of the 
planes Zi into the trace formula; the contributions of 
longer orbits were shown to be negligible Q . 

In the upper parts of Figs. |l| and ||, the shape bound- 
aries p{z) are shown for two deformations with fixed elon- 
gation c = 1.53. The first corresponds to the symmetric 
outer barrier with a = 0, and the second to the asymmet- 
ric saddle with a = 0.13. The vertical lines show the po- 
sitions Zi of the planes containing the periodic orbits. In 
the symmetric case (Fig. 1^) , the shortest orbits at zq = 
are slightly unstable, wheras those at —Z2 — z\ = 0.414 
are stable. In the asymmetric case (Fig. H), there is only 
one plane at zq = 0.671 containing stable orbits. In the 
lower parts of the figures, we show Poincare surfaces of 
section {p, (jj) for trajectories with (conserved) angular 
momentum component L^ = 0, where p is the compo- 
nent of the momentum parallel to the tangent plane at 
a reflection point and (p is the polar angle in the {z,p) 
plane of a reflection point at the boundary. (Note that 
the Poincare mapping in these variables is not area pre- 
serving; this is, however, immaterial for the qualitative 
interpretation of the figures.) In both cases we find two 
main islands of stability containing the fixed points corre- 
sponding to the shortest stable and unstable orbits men- 
tioned above. Apart from some KAM chains of smaller 
islands corresponding to higher resonances, the remain- 
der part of the phase space is mainly chaotic. Note that 
the overall chaoticity is much larger in the energetically 
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more stable asymmetric case. We thus find that the shell 
effect which leads to the energetic instability of the sym- 
metric outer barrier, at the same time pulls the system 
into a more chaotic transition state. 

In Ref. PI, the microscopic origin of the instability 
against asymmetric deformations was linked to specific 
quantum states having their probability maxima in the 
central equatorial plane or in two planes parallel to it. 
These states were energetically most sensitive to the 
asymmetric deformations, whereas all others were prac- 
tically not. In order to check this behaviour in our 
present cavity model and to investigate the relation of 
these states to the leading periodic orbits, we have cal- 
culated its quantum spectrum (with Dirichlet boundary 
conditions) at the relevant deformations. In Fig. y we 
show the wave number spectrum ki = \/2mEi/h (in units 
of l/R) for the lowest states i with angular momentum 
Lz — 4: obtained at c = 1.53 versus the asymmetry pa- 
rameter a. Whereas most levels vary only little with a, 
we clearly recognize a sequence of states whose energies 
decrease more and more steeply with increasing number 
i. These states are diabatic in that they connect different 
portions of the level spectrum through avoided crossings; 
in the upper part of the spectrum we have emphasized 
them by heavy grey lines and labeled them by the num- 
bers i of the adiabatic levels counted from the bottom. 
We shall show below that the energies of these diabatic 
states can be obtained from an EBK quantization of the 
classical motion near the equatorial planes. 

The mechanism of the energy instability against a now 
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FIG. 2. The same as in Fig. hi for asymmetry a = 0.13. 



3. Quantum levels kiR for the lowest states i with 
Lz — 4 ai elongation c — 1.53, plotted versus asymmetry 
a. The emphasized portions connected by avoided crossings 
denote some of the diabatic states (see text). 



is this : if one or two of the diabatic states are located 
just underneath the Fermi energy (and thus occupied), 
the system gains energy when a is increased from zero. 
Since all the occupied "inert" states slightly increase at 
larger values of a, the sum of occupied levels (and hence 
5E, see Ref. M) will exhibit a minimum at some finite 
value Qfo- In our model calculation MM for Pu^^° with a 
fixed Fermi energy kp = 12.1/i?, we obtained ao ~ 0.13. 

In Fig. Hwe show probability distributions by contour 
plots of the squares of wave functions \tpi{z,p)\'^ in the 
{z,p) plane. The centre and bottom panels correspond 
to i = 10 at a = and to i — 7 &i a — 0.2, respectively, 
representing one diabatic state at the two ends of the a 
interval. The top panel represents the beginning of the 
next higher diabatic state, starting asi = llatQ; = 0. 
The probability maxima are clearly located precisely in 
the planes Zi of the shortest periodic orbits, indicated by 
the heavy vertical lines. This pattern was found to be 
consistent: the probability maxima of all diabatic states 
are located in the planes of the shortest periodic orbits, 
whereas the maxima of all other "inert" states could not 
be correlated to any of the leading periodic orbits. 

In order to quantify the correspondence between the 
diabatic quantum states and the leading periodic orbits, 
we apply a method analogous to the quasimode construe- 
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FIG. 4. Contour plots of the probability distributions 
\ipi{z,p)'^ for selected states with L^ = 4 at the elongation 
c = 1.53. The heavy black lines show the shape p[z) of the 
boundary. The grey vertical lines indicate the planes Zi con- 
taining the shortest periodic orbits. Top: state i = 11 at 
asymmetry a = 0.0. Centre: state i = 10 at asymmetry 
a = 0.0. Bottom: state i = 7 at asymmetry a = 0.2. 



tion for eigenfunctions located in the vicinity of stable 
periodic orbits (see, e.g., fl^l)- Consider the classical mo- 
tion in the vicinity of a stable equatorial plane. It consists 
of small oscillations around this plane. In a linearized ap- 
proximation, the motion is integrable and restricted to 
tori in phase space. The quasimode approximation con- 
sists in applying the EBK quantization conditions to this 
torus structure pl| . To derive the quantization condi- 
tions, we use the fact that the quasimode approximation 
depends only on the two curvature radii at the boundary 
of the stable orbit plane: the radius Ri in the equatorial 
plane itself, and the radius i?2 perpendicular to it. We 
therefore replace our billiard system locally by an axially 
symmetric ellipsoid with the same curvature radii of the 
equatorial planes. For an ellipsoid, the classical action 
variables are known [[l2| . In the vicinity of the equatorial 
plane they are given by the z component of the angular 
momentum, Lz/ti — 0, ±1, ±2, . . . , by the action Ip for 
the radial motion in the equatorial plane. 
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p^Rl ^ LI - \Lz\aTCCos\Lz/ipRi)\ 



K arcsin yjRi/R2 — L'^ / {p'^ R1R2) 
2 ^pmiiR2-Ri)+Ll 



(1) 



and by the action I^ for the motion out of this plane, 

/. « n/[2^p^RiiR2-Ri) + Ll] . (2) 

These actions depend on the three constants of motion: 
the energy E (or the momentum squared p^ = % /c^), 
the angular momentum L^, and a conserved quantity n 
which vanishes for the motion in the equatorial plane 
and is the analogue of the second constant of motion in 
a two-dimensional elliptic billiard (given by the product 
of angular momenta with respect to the two foci). The 
actions (1^), (0) are already given here in the linearized 
approximation obtained from the exact actions (given by 
elliptic integrals) through Taylor expansion up to first 
order in n. Eliminating k from Eqs. (pj), (0), we arrive at 



Ip = 



p^Rf - LI - \L,\arccos\Lz/{pRi)\ 



h aTCsmy/Ri/R2 - Ll/{p^RiR2) , (3) 



to be used with the EBK quantization conditions 



h{n- 1/2), 
h (m - 1/4) . 



n = 

TTL 



1,2,3,. 
= 1,2,3,. 



(4) 
(5) 



Note that if one puts Iz — above, Eqs. M) and (ra) yield 
precisely the radial EBK quantization condition (without 
further approximation) for a two-dimensional circular bil- 
liard with radius Ri. 

Equations (0) - (m represent our quasimode approxi- 
mation. They implicitly yield the energies Enm (or wave 



numbers fc„m) of quantum states with angular momen- 
tum Lz located near the planes of stable periodic orbits, 
which we expect to represent the diabatic states discussed 
above. The quantum numbers n and m, count the wave 
function extrema in the z and p directions, respectively. 
(Correspondingly, n—\ and m— 1 count the numbers of 
nodes in the respective directions.) E.g., the squared 
wavefunctions shown in the centre and bottom panels of 
Fig. 1^ should represent the state (1,2), and the one shown 
in the top panel should belong to the state (2, 2). 

As a test of our interpretation of the diabatic quan- 
tum states and their semiclassical quantization, we com- 
pare in Fig. g the exact quantum spectrum hi (shown 
by solid lines) to the approximate EBK levels knrm cal- 
culated here for the lowest states with angular momen- 
tum Lz = 0. Due to our linearization of the actions (|^), 
(0) , the agreement should be best for small values of the 
quantum number n. It is, indeed, perfect for all states 
with n = 1 shown in Fig. H by the dotted lines, which 
agree exactly with the quantum levels of the correspond- 
ing diabatic states. The agreement is less good for the 
states with n = 2 (shown by the dashed lines), but our 
semiclassical approximation still reproduces their correct 
slopes for larger a and allows for their unique assignment. 
Similar results are also obtained for larger values of Lz 
(for which the numerical agreement actually improves). 

In summary, we have established a correspondence be- 
tween the shortest classical periodic orbits and a set of 
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FIG. 5. Quantum levels kiR of lowest states i with L^ — 
at elongation c = 1.53 versus asymmetry a (solid lines). La- 
bels {n,m) indicate diabatic states; their levels k„mR ob- 
tained by the quasimode approximation, Eqs. (|3|) - (n), are 
shown by dotted lines for n = 1 and dashed lines for n = 2. 



diabatic quantum states responsible for the onset of the 
left-right asymmetry in a simple axially deformed cav- 
ity model describing schematically the shapes of an ac- 
tinide nucleus on its adiabatic path to fission. The energy 
gain due to the asymmetric deformations near the outer 
barrier had earlier been reproduced by a semiclassical 
trace formula [b|J4| ; thereby the fission path through the 
deformation space was shown to be determined by the 
constancy of the actions of the shortest periodic orbits. 
In the present paper we have demonstrated that these 
orbits lie at the centres of small islands of stability in 
an otherwise chaotic phase space. Its degree of chaotic- 
ity is increased by the shell effect causing the asymmet- 
ric deformations. We have also shown that the diabatic 
quantum states which energetically favor the asymme- 
try have their probability maxima precisely in the planes 
where the shortest periodic orbits are located. A quasi- 
mode approximation based on EBK quantization of the 
linearized classical motion around the planes of the short- 
est stable orbits allowed us to uniquely assign quantum 
numbers {Lz,n,m) to the diabatic quantum states and 
to semiclassically reproduce their energies rather well. 

We acknowledge the help of P. Meier in constructing 
the energy level plots. This work has been supported by 
the Deutsche Forschungsgemeinschaft (M.B.) and a Ha- 
bilitationsstipendium des Freistaates Bayern (S.M.R.). 
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